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Abstract 

In  this  paper,  we  develop  a  reduced  basis  construction  method  that  allows  for  separate  consideration  of  baseline  and  actuated 
dynamics  in  the  reduced  modeling  process.  A  prototype  initial  boundary  value  problem,  governed  by  the  two-dimensional 
Burgers  equation,  is  formulated  to  demonstrate  the  utility  of  the  method  in  a  boundary  control  setting.  Comparisons  are 
done  between  reduced  and  full  order  solutions  under  open-loop  boundary  actuation  to  illustrate  advantages  gained  by  separate 
consideration  of  actuated  dynamics.  A  tracking  control  problem  is  specified  using  a  linear  quadratic  regulator  formulation. 
Comparisons  of  feedback  control  effectiveness  are  done  to  demonstrate  benefits  in  control  effectiveness  obtained  from  separate 
consideration  of  actuated  dynamics  during  model  reduction. 


Nomenclature 


L 

S 

N 

T 

Na 
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Mb  ,  Ma 
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A,  v 
A,  B 
G,F 
Q,R 
7 

n 

t 
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covariance  matrix 

solution  snapshot 

total  number  of  snapshots 

actuated  snapshot 

total  number  of  actuated  snapshots 

basis  mode  and  temporal  coefficient,  respectively 

number  of  basis  modes 

number  of  baseline  and  actuator  modes,  respectively 

baseline  and  actuator  mode,  respectively 

baseline  and  actuator  mode  temporal  coefficient,  respectively 

eigenvalue  and  eigenvector,  respectively 

state  and  control  matrices,  respectively 

nonlinear  and  forcing  matrices,  respectively 

state  and  control  weight  matrices,  respectively 

control  robustness  parameter 

spatial  domain 

time 

spatial  coordinate  vector 
spatial  step-size 


Introduction 


Reduced  order  modeling  has  received  significant  research  attention  in  recent  years.  For  many  problems 
of  practical  interest,  the  order  of  the  system  describing  the  application  must  be  reduced.  An  illustrative 
example  where  this  is  required  is  the  development  of  feedback  control  laws  for  fluid  flow  configurations.  It 
is  not  uncommon  for  discretized  flow  models  to  describe  millions  of  state  variables.1  Unfortunately,  the 
development  of  systematic  feedback  control  laws  from  systems  of  such  large  dimension  is  a  computationally 
intractable  problem.  For  example,  if  one  uses  a  linear  quadratic  regulator  (LQR)  control  formulation,  roughly 
1012  Riccati  unknowns  need  to  be  calculated  for  a  discretized  flow  model  describing  106  states.  The  Riccati 

1This  paper  is  declared  a  work  of  the  U.S.  Government  and  is  not  subject  to  copyright  protection  in  the  United  States. 
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unknowns  are  solutions  to  a  nonlinear  matrix  equation.2  Existing  computing  power  and  computational 
algorithms  are  not  capable  of  solving  an  LQR  problem  of  such  large  dimension.  For  dynamical  models  that 
are  very  large  scale,  such  as  those  describing  fluid  flow  configurations,  it  is  apparent  that  the  order  of  the 
system  must  be  reduced  prior  to  control  law  design.3 

Several  research  efforts  have  been  concerned  with,  or  relied  upon,  the  development  of  order  reduction 
strategies  that  provide  reduced  order  models  in  a  form  amenable  to  state-space  feedback  control  law  de¬ 
sign. 4-14  For  the  case  of  boundary  control,  the  development  of  such  models  has  been  an  open  problem. 
For  many  applications  of  practical  interest,  such  as  feedback  control  of  the  air  flow  over  an  airplane  wing, 
boundary  actuation  is  a  requirement.  These  applications  require  actuation  to  be  located  on  the  surface  if 
they  are  to  be  implemented  in  hardware  in  the  physical  system.  Most  reduced  modeling  efforts  have  either 
been  concerned  with  the  case  of  control  action  via  a  body  force  or  have  approximated  boundary  actuation 
by  a  body  force  on  the  domain  interior.  Few  techniques  have  been  obtained  thus  far  providing  reduced  order 
models  with  explicit  boundary  input  that  can  be  used  for  systematic  design  of  feedback  control  laws.  It  is 
difficult  to  extract  the  boundary  control  input  in  the  reduced  model.  When  proper  orthogonal  decomposition 
(POD)  is  used  in  conjunction  with  Galerkin  projection,  boundary  conditions  are  absorbed  in  the  process. 
An  additional  step  is  required  to  make  the  action  of  the  control  input  explicit  in  the  model. 

Recent  inroads  have  allowed  for  accurate  implementation  of  boundary  control  in  reduced  order  state-space 
models.  One  approach  is  to  capture  the  influence  of  boundary  actuation  indirectly.  Of  note  is  the  work  of 
Serrani  and  Kasnakoglu15  where  a  basis  representing  the  influence  of  boundary  actuation  is  captured  by  solv¬ 
ing  a  particular  optimization  problem  in  L2.  The  dimensions  of  reduced  order  dynamical  models  constructed 
with  this  method  are  typically  very  small,  allowing  for  implementation  of  nonlinear  control  techniques  where 
small  system  order  is  critical.  Other  methods  utilize  POD  in  combination  with  a  weak  formulation  of  the 
Galerkin  projection. 16-19  The  advantage  of  a  weak  formulation  is  that  boundary  conditions  appear  naturally 
when  the  reduced  model  is  written  weakly.  Posing  the  system  weakly  also  reduces  regularity  requirements 
of  the  solution,  reducing  errors  that  result  from  numerical  approximation  of  high  order  derivatives  in  the 
Galerkin  system.  In  this  paper,  we  utilize  a  method16  incorporating  difference  approximations  in  the  weak 
formulation  to  construct  reduced  models  with  boundary  control  input  appearing  explicitly. 

Finally,  order  reduction  methods  must  be  flexible  so  that  they  are  amenable  to  either  simulation  or 
experimental  approaches.  Each  approach  has  advantages  and  disadvantages,  but  the  particular  application 
at  hand  often  dictates  the  approach  taken.  Methods  solely  suitable  for  a  simulation  approach  will  not  likely 
lend  themselves  to  implementation  in  an  experimental  facility.  There  is  currently  a  need  for  order  reduction 
strategies,  suitable  for  practical  control  configurations,  that  can  utilize  system  data  generated  by  simulation 
or  experiment.  It  is  in  this  spirit  that  the  techniques  presented  in  this  paper  are  developed.  Simulation  data 
is  used  to  demonstrate  methods  in  the  current  paper.  The  application  of  these  methods  to  an  experimental 
laser  turret  configuration  is  presented  in  a  companion  paper.20 

Proper  Orthogonal  Decomposition 


Proper  orthogonal  decomposition  is  a  popular  technique  used  to  construct  optimal  basis  functions  attrac¬ 
tive  for  reduced  order  modeling.  In  this  work,  we  use  a  POD  algorithm21  based  on  the  snapshot  method22 
to  construct  the  low  order  basis  needed  for  the  development  of  the  reduced  order  model.  For  the  sake  of 
completeness  and  clarity  in  what  follows,  we  briefly  describe  this  technique.  A  data  ensemble  of  snapshots 
{S',(x)}^.1  is  generated  for  the  system  via  numerical  simulation  or  experiment,  where  N  is  the  total  number 
of  snapshots.  Each  snapshot  consists  of  instantaneous  system  data.  With  the  snapshot  ensemble  in  hand, 
the  N  x  N  correlation  matrix  L  defined  by 


Li,j  =  (Si,  Sj) 

is  constructed.  In  this  work,  we  utilize  the  standard  L2(fl)  inner  product 

(Si,Sj)=  [  SiS* dx, 

Jn 

where  S*  denotes  the  complex  conjugate  of  Sj,  in  the  construction  of  L. 


(1) 

(2) 
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The  eigenvalues  {Aj}£Lj  of  L  are  calculated  and  sorted  in  descending  order.  The  ratio 

a,  \ 


100 


Xf=i^ 


(3) 


is  used  to  determine  the  number  of  POD  basis  functions  to  construct.  The  quantity  in  (3)  provides  a  measure 
of  the  ensemble  energy  that  is  captured  by  a  POD  basis  consisting  of  M  modes.  By  requiring  a  percentage 
of  the  energy  contained  in  the  snapshot  ensemble  be  contained  in  the  basis,  the  smallest  value  of  M  is 
calculated  such  that  the  quantity  in  (3)  is  greater  than  or  equal  to  that  percentage. 

The  eigenvectors  {uj}i=1  corresponding  to  the  M  eigenvalues  of  largest  magnitude  are  calculated.  Each 
eigenvector  is  normalized  so  that 


(4) 


The  orthonormal  POD  basis  set  {^(x)}^1  is  constructed  according  to 


N 

<MX)  = 

i= i 


(5) 


where  vtj  is  the  jth  component  of  Uj.  With  the  basis  in  hand,  the  system  solution  w[t,  x)  is  approximated 
as  a  linear  combination  of  POD  modes,  i.e., 


M 

w(t,x)  rts  (6) 

i— 1 

Split-POD 

Using  an  energy  argument  based  on  (3)  is  a  convenient  way  to  determine  the  number  of  POD  modes 
needed  for  the  reduced  model.  However,  blindly  applying  (3)  to  determine  the  necessary  number  of  modes 
is  problematic  in  a  boundary  control  setting.  In  boundary  control  applications,  it  is  desired  that  control 
input  energy  be  as  small  as  possible  while  still  satisfying  the  control  objective.  For  example,  in  flow  control 
applications  where  the  control  is  located  on  a  surface,  it  is  desired  that  small  control  inputs  yield  large 
changes  in  the  flow  field  behavior.23  In  essence,  effects  of  small  control  inputs  are  amplified  by  the  natural 
instabilities  and  comparatively  high  energy  content  of  the  baseline  flow  field.  Simply  applying  (3)  to  an 
ensemble  consisting  of  baseline  and  actuated  data  presents  the  risk  of  important  structures  due  to  control 
input  being  discarded.  These  structures  are  typically  of  much  lower  energy  content  than  those  associated 
with  the  baseline  solution.  For  these  reasons,  we  extend  the  POD  algorithm  described  above  so  that  baseline 
and  control  input  energy  are  considered  separately.  We  refer  to  this  method  as  split-POD.  The  basis  resulting 
from  this  method  will  consist  of  modes  significant  to  the  baseline  solution  as  well  as  those  significant  from 
an  actuation  standpoint. 

The  basic  idea  is  to  decompose  each  snapshot  in  the  ensemble  into  a  component  in  the  span  of  a  baseline 
POD  basis  and  an  orthogonal  component.  This  is  done  by  employing  useful  properties  of  orthogonal  projec¬ 
tions  on  Hilbert  spaces.24  By  considering  the  case  of  baseline  and  actuated  data  separately,  the  orthogonal 
component  is  constructed  such  that  it  contains  new  information  due  to  the  control  input. 

An  ensemble  of  solution  snapshots  is  generated  for  the  case  of  zero  control  input.  From  this  baseline 
snapshot  ensemble,  a  set  of  POD  modes  is  constructed  that  contains  a  high  percentage  of  the  baseline  energy. 
For  notational  convenience,  denote  this  baseline  basis  by  {£,■ where  Mb  is  the  number  of  modes.  By 
employing  the  energy  ratio  in  (3),  we  choose  Mb  so  that  the  baseline  basis  contains  most  of  the  energy 
contained  in  the  baseline  snapshot  ensemble. 

With  a  large  set  of  baseline  modes  in  hand,  an  ensemble  of  solution  snapshots  is  generated  for  the  case 
of  nonzero  control  input.  Denote  the  actuated  ensemble  by  {Tj}^  where  N a  is  the  number  of  actuated 
snapshots. 

For  each  snapshot  in  the  actuated  ensemble,  we  determine  the  component  that  is  in  the  span  of  the 
baseline  basis.  In  particular,  define  according  to 

bij  =  (Ti,  &>,!<*<  Na,  1  <j<MB-  (7) 


3 


Then,  bij  is  the  projection  of  the  z'tli  actuated  snapshot  T.t  onto  the  jth  baseline  POD  mode  £j.  In  other 
words,  the  product  bij£j  is  the  component  of  T,;  that  is  in  the  direction  of  £j.  The  linear  combination 
bij£j  is  the  component  of  T,  in  the  span  of  the  baseline  basis. 

Define  Tt  according  to 

Mb 

Ti  =  Ti  —  'y  ~  bij^j.  (8) 

3= i 

Then,  Tj  is  the  component  of  T,;  not  contained  in  the  span  of  the  high  order  baseline  basis.  As  T,;  is  a  solution 
snapshot  for  the  case  of  nonzero  control  input,  T,  consists  of  new  information  due  to  the  control  input.  A 
second  set  of  POD  modes  is  constructed  from  the  data  ensemble  Denote  this  set  of  “actuator 

modes”  by  where  M a  is  the  number  of  modes.  The  energy  ratio  in  (3)  is  used  to  determine  Ma 

such  that  the  basis  of  actuator  modes  contains  an  arbitrary  amount  of  the  additional  energy  resulting  from 
the  control  input.  We  have  employed  the  high  order  baseline  basis  to  extricate  new  information  resulting 
from  the  control.  Having  done  that,  the  baseline  basis  is  now  truncated  to  contain  a  more  reasonable  number 
of  modes  for  model  development.  Denote  Mb  as  the  number  of  baseline  modes  retained  after  truncation 
with  corresponding  basis  set  {£,■  }j2i 

To  simplify  reduced  order  modeling  via  Galerkin  projection,  it  is  advantageous  to  use  a  basis  consisting 
of  orthonormal  modes.  We  now  show  that  the  sets  of  baseline  and  actuator  modes  can  be  combined  into  an 
overall  basis  set  where  all  modes  are  orthonormal. 

By  construction,  baseline  modes  are  orthonormal.  Similarly,  actuator  modes  are  orthonormal.  Consider 
the  inner  product  of  and  Tt  for  arbitrary  i,  k.  We  see  that 

/mb  \ 

(Ti,  £fc)  =  \  ~  'y  'j  bij£j,  £,kj 

Mb 

=  (Ti,  £fc)  —  b^  (£j,  £&} 

3= i 

As  the  baseline  modes  are  orthonormal,  {£j,f;k)  =  0  unless  j  =  k.  For  j  =  k,  (£j,£k)  =  1-  As  a  result, 

mb 

{Ti,£k)  -  y {€j>€k)  =  {Ti,£k)  ~  bik-  (11) 

3=1 

By  definition  (7),  bn-  =  (Tj,^).  Thus,  we  have 

=  0,  for  arbitrary  i,k.  (12) 

From  (5),  each  actuator  mode  rji  is  a  linear  combination  of  the  snapshots  {Tj}^.  As  a  result,  (12)  in 
combination  with  the  linearity  of  the  inner  product  yields 

(Vit^k)  —  0,  for  arbitrary  i,k.  (13) 

Therefore,  we  have 

{Vi}  A  {£,}•  (14) 

This  result  allows  us  to  combine  the  baseline  and  actuator  modes  into  an  overall  basis  set 

(15) 

where  all  modes  in  the  basis  are  orthonormal.  The  system  solution  w{t,  x)  is  still  approximated  as  a  linear 
combination  of  modes  as  in  (6).  Moreover,  separate  consideration  of  baseline  and  actuated  energy  allows  us 
to  write  this  linear  combination  as 

Mb  Ma 

w(t,x)  «  5^  0J  (*)£,- (x)  +5^A(t)»7i(x),  (16) 

3=1  *= 1 


(9) 

(10) 
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where  { 0:] }^f_”  and  {A}^  are  temporal  coefficients  for  the  baseline  and  actuator  basis,  respectively.  This 
allows  us  to  consider  the  system  solution  as  a  baseline  component  and  an  additional  component  induced  by 
the  control  input. 

Model  Problem 


We  now  demonstrate  the  benefit  of  careful  consideration  of  actuated  dynamics  on  reduced  model  accuracy 
and  subsequent  boundary  control  effectiveness.  A  distributed  parameter  system  is  formulated  which  models 
convective  flow  over  an  obstacle.  Let  Q-|  C  M2  be  the  rectangle  given  by  (a,  b]  x  ( c,d ).  Let  Q2  C  17 1  be  the 
rectangle  given  by  [ai,  02]  x  [61, 62]  where  a  <  a±  <  a2  <  b  and  c  <  61  <  62  <  d.  The  problem  domain,  17,  is 
given  by  17  =  fli  \  172-  In  this  configuration,  H2  is  the  obstacle.  Dirichlet  boundary  controls  are  located  on 
the  obstacle  bottom  and  top,  denoted  by  Ts  and  Tt,  respectively. 


Figure  1:  Problem  Geometry. 


The  dynamics  of  the  system  are  described  by  the  two-dimensional  Burgers  equation 


d  1 

w(t,  x,  y)  +  V  •  F(w)  =  —Aw(t,  x,  y) 


dt 

for  t  >  0  and  (x,y)  £  17.  In  (17),  F(w)  has  the  form 

w2(t,x,y) 


Re 


F(w)  = 


Ci- 


C' 


w2{t1  x,  y) 


1  T 


(17) 


(18) 


where  C\,  C2  are  nonnegative  constants.  This  equation  has  a  convective  nonlinearity  like  that  found  in 
the  Navier-Stokes  momentum  equation  modeling  fluid  flow.25  The  quantity  Re,  a  nonnegative  constant,  is 
analogous  to  the  Reynolds  number  in  the  Navier-Stokes  momentum  equation. 

To  complete  the  model  of  the  system,  boundary  conditions  must  be  specified  as  well  as  an  initial  condition. 
For  simplicity,  boundary  controls  are  assumed  to  be  separable.  With  this  assumption,  we  specify  conditions 
on  the  obstacle  bottom  and  top  of  the  form 


w(t,TB)  =  uB(t)^B(x),  (19) 

w(t,VT)  =  uT(t)^fT(x).  (20) 

In  (19)- (20),  uB(t)  and  ur{t)  are  the  controls  on  the  bottom  and  top  of  the  obstacle,  respectively.  The  profile 
functions  $5(1)  and  T?^#)  describe  the  spatial  influence  of  the  controls  on  the  boundary.  A  parabolic  inflow 
condition  is  specified  of  the  form 

w(t,Tin)  =  f(y).  (21) 

At  the  outflow,  a  Neumann  condition  is  specified  according  to 

^w(t,Tout)  =  0.  (22) 

For  notational  convenience,  denote  the  remaining  boundary  as  Tj/.  We  require  that  values  be  fixed  at  zero 
along  Tjj  as  time  evolves.  The  resulting  boundary  condition  is  of  the  form 

«j(t,r^)  =  0.  (23) 
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The  initial  condition  of  the  system  is  given  by 

w(0,x,y)  =  w0(x,y)  £  L2(n).  (24) 

As  developed  in  previous  work,1'  the  reduced  order  state-space  model  for  this  system  with  explicit 
boundary  control  input  is  of  the  form 

a  =  Aa  +  Bu+  G(a)  +  F.  (25) 

Projecting  the  initial  condition  Wo(x,y)  onto  the  POD  basis  results  in  an  initial  condition  for  the  reduced 
order  model  of  the  form 

a(0)  =  a0.  (26) 

Open-Loop  Comparisons 

We  demonstrate  the  impact  of  basis  construction  technique  on  the  ability  of  the  reduced  model  to 
represent  dynamics  induced  by  a  boundary  control  input.  Instantaneous  snapshots  are  generated  for  (17), 
(19)-(24)  via  numerical  simulation.  A  positive  parabolic  profile  with  unit  maximum  amplitude  is  specified 
for  the  inlet  condition  in  boundary  condition  (21).  In  (18),  we  set  C i  =  1  and  =  0  in  order  to  obtain 
solutions  that  convect  from  left  to  right  for  the  positive  inlet.  In  addition,  we  specify  that  Re  =  300.  The 
problem  domain  O  is  discretized,  resulting  in  a  uniform  grid  with  spatial  step-size  h.  We  utilize  a  finite 
difference  scheme26  to  numerically  solve  the  model  problem  with  and  without  boundary  control  input.  The 
resulting  discretized  system  describes  roughly  2,000  states. 

We  construct  reduced  order  models  from  snapshot  ensembles  obtained  for  two  scenarios.  In  the  first 
scenario,  snapshots  are  generated  for  the  baseline  solution  and  for  boundary  actuation  at  a  fixed  frequency. 
For  the  second  case,  snapshots  are  generated  for  the  baseline  solution  and  a  more  complicated  chirp  input 
where  the  frequency  varies  with  time.  For  both  scenarios,  we  compare  model  agreement  resulting  from 
combining  baseline  and  actuated  snapshots  into  an  overall  lumped  snapshot  ensemble  to  that  obtained  by 
decomposing  snapshots  into  their  baseline  and  actuated  components  and  constructing  the  split-POD  basis 
as  in  (15).  In  the  results  that  follow,  basis  construction  from  a  lumped  snapshot  ensemble  is  referred  to 
as  the  “lumped”  method.  Constructing  the  basis  by  considering  baseline  and  actuated  energy  separately  is 
referred  to  as  the  “split”  method. 

Scenario  1 

Snapshots  are  generated  for  the  baseline  solution  and  for  the  solution  arising  under  periodic  boundary 
actuation.  Inputs  specified  are  of  the  form 

Us{t)  =  sin(nt)  Ur(t)  —  0,  (27) 

us{t)  =  0  UT{t)  =  sin(irt).  (28) 

In  (27),  periodic  actuation  is  done  on  the  bottom  of  the  obstacle  while  values  along  the  obstacle  top  are  held 
at  zero.  In  (28),  values  along  the  obstacle  bottom  are  held  at  zero  with  periodic  actuation  occurring  at  the 
top.  For  each  control  input  listed  in  (27)-(28),  snapshots  are  taken  in  increments  of  At  =  0.1  starting  from 
t  =  0  and  ending  at  f  =  15.  The  steady  baseline  solution  is  used  for  the  initial  condition. 

With  an  ensemble  of  snapshots  in  hand,  ratio  (3)  is  used  to  determine  the  number  of  basis  modes  to 
construct.  Requiring  that  99.9%  of  the  ensemble  energy  be  contained  in  the  POD  basis  results  in  a  lumped 
basis  consisting  of  7  modes.  Separate  consideration  of  baseline  and  actuated  energy  results  in  a  split  basis 
consisting  of  1  baseline  mode  and  16  actuator  modes. 

We  now  employ  linear  combination  (6)  to  compare  boundary  condition  agreement  between  the  full  order 
system  and  reduced  models  obtained  via  the  lumped  and  split-POD  bases.  By  specifying  characteristic 
functions  for  the  control  profiles  $b(j;)  and  T^a;)  in  (19)-(20),  we  see  that 

M 

B)  «  io(t,  Tb)  =uB{t),  (29) 

i= 1 

M 

5>(*)*(Ft)  «  w{t,  rT)  =  uT{t),  (30) 

i= 1 
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We  construct  the  linear  combinations  on  the  left  in  (29)- (30)  and  compare  them  to  the  exact  boundary 
conditions  us{t )  and  ur(t)  specified  in  the  full  order  system.  We  first  perform  comparisons  for  boundary 
inputs  explicitly  used  during  ensemble  generation.  Results  obtained  for  the  baseline  solution  and  for  the 
solution  with  periodic  boundary  actuation  of  the  form  sm(-7rf)  are  shown  in  Figure  2.  In  that  figure,  dashed 


Figure  2:  Boundary  Condition  Accuracy. 


curves  denote  the  linear  combination  of  POD  modes  restricted  to  the  boundary.  Solid  curves  denote  the 
exact  full  order  boundary  input.  Results  obtained  for  the  lumped  basis  method  are  plotted  on  the  left. 
Split-POD  basis  results  are  plotted  on  the  right.  As  seen  in  Figure  2,  both  methods  result  in  very  good 
agreement  between  the  exact  boundary  conditions  and  the  linear  combination  of  POD  modes  restricted  to 
the  boundary.  Dashed  and  solid  curves  are  virtually  identical. 

In  a  feedback  control  setting,  dynamics  induced  by  the  control  are  typically  not  known  a  priori.  The 
specifics  of  the  control  law  and  the  dynamics  induced  by  it  are  not  known  prior  to  ensemble  creation. 
Typically,  in  the  closed-loop  system,  the  boundary  input  resulting  from  the  control  law  will  be  different 
than  the  inputs  used  to  generate  the  snapshot  ensemble.  As  a  result,  it  is  useful  to  compare  reduced  and 
full  order  model  agreement  for  inputs  not  specified  during  ensemble  creation.  This  provides  insight  into 
the  suitability  of  the  reduced  model  for  closed-loop  control  law  design.  For  these  reasons,  we  now  compare 
boundary  condition  agreement  between  the  reduced  and  full  order  systems  for  open-loop  inputs  that  were 
not  used  during  ensemble  creation.  Boundary  inputs  specified  are  of  the  form 


usit)  =  min 


uT{t) 


(31) 

(32) 


Results  obtained  for  the  lumped  and  split  basis  methods  are  shown  in  Figure  3. 

As  seen  in  Figure  3,  very  good  agreement  is  seen  between  the  reduced  and  full  order  systems  for  the 
split  method,  even  though  the  inputs  considered  were  not  specifically  included  in  the  snapshot  ensemble. 
Condition  (32)  is  reconstructed  well  using  the  lumped  method.  However,  the  reconstruction  of  the  piecewise 
linear  input  in  (31)  is  much  less  accurate  when  the  lumped  basis  is  used. 

To  further  compare  the  lumped  and  split  basis  methods  and  their  utility  for  control  law  design,  we  project 
the  full  order  solution  at  each  time  step  onto  the  lumped  and  split  POD  bases.  The  resulting  temporal 
coefficients  are  compared  to  those  predicted  by  the  reduced  order  models.  Boundary  inputs  specified  are 
as  in  (31)-(32).  Results  obtained  for  the  first  5  temporal  coefficients  of  the  lumped  method  are  shown  in 
Figure  4.  The  first  5  temporal  coefficients  for  the  split  method  are  shown  in  Figure  5.  In  Figures  4-5, 
solid  curves  denote  values  of  temporal  coefficients  obtained  via  the  projection.  Dashed  curves  denote  the 
solution  of  the  reduced  order  model.  Overall,  the  accuracy  of  the  split  method  is  better,  particularly  for 
temporal  coefficients  with  higher  frequency  content.  Separate  consideration  of  actuated  energy  during  basis 
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Lumped  Split 


Figure  3:  Boundary  Condition  Accuracy. 


Lumped 


Figure  4:  Temporal  Coefficient  Accuracy  for  the  Lumped  Method. 


construction  results  in  better  representation  of  dynamics  induced  by  boundary  inputs  not  specified  during 
ensemble  creation.  This  is  advantageous  in  a  boundary  feedback  control  setting  where  dynamics  induced  by 
the  control  are  typically  not  known  beforehand. 

Scenario  2 


It  is  likely  that  a  snapshot  ensemble  for  inputs  at  a  single  frequency  results  in  a  POD  basis  that  does 
not  adequately  span  the  dynamics  induced  by  a  feedback  control.  Feedback  controls  designed  from  such  a 
basis  are  bound  to  be  ineffective  when  implemented  in  the  full  order  system.  In  Scenario  2,  we  compare  the 
lumped  and  split  basis  methods  using  a  snapshot  ensemble  generated  from  chirp  inputs  of  the  form 


us{t)  =  sin  (w  (e*)0 

uT(t)  =  0, 

(33) 

UB(t)  =  0 

ut{S)  =  sin  (^n  (e*)°'3^  . 

(34) 

As  seen  in  Figure  6,  an  input  of  this  form  generates  the  system  response  over  a  range  of  frequencies.  The 
resulting  POD  basis  is  much  more  likely  to  sufficiently  span  the  unknown  dynamics  generated  by  a  feedback 
control. 


Split 


Figure  5:  Temporal  Coefficient  Accuracy  for  the  Split  Method. 


Instantaneous  snapshots  are  generated  for  the  baseline  solution  as  well  as  for  solutions  arising  from  the 
inputs  in  (33)- (34).  For  each  control  input  listed  in  (33)- (34),  snapshots  are  taken  in  increments  of  At  =  0.1 
starting  from  t  =  0  and  ending  at  t  =  15.  The  steady  baseline  solution  is  used  for  the  initial  condition. 

Requiring  that  99.9%  of  the  ensemble  energy  be  contained  in  the  POD  basis  results  in  a  lumped  basis 
consisting  of  7  modes.  A  split  basis  comprised  of  1  baseline  mode  and  25  actuator  modes  are  needed  when 
baseline  and  actuated  energy  is  considered  separately. 

As  in  Scenario  1,  we  compare  boundary  condition  accuracy  of  the  lumped  and  split  methods  for  boundary 
inputs  not  specified  during  ensemble  generation.  For  the  sake  of  comparison,  we  use  the  boundary  conditions 
given  by  (31)-(32).  Results  obtained  for  the  lumped  and  split  basis  methods  are  shown  in  Figure  7.  As  seen 
in  that  figure,  very  good  boundary  condition  agreement  is  seen  between  the  reduced  and  full  order  systems 
for  both  basis  methods.  In  particular,  by  comparing  Figures  3  and  7,  we  see  that  the  piecewise  linear 
boundary  condition  in  (31)  is  represented  much  better  by  the  lumped  method  when  inputs  (33)- (34)  are 
used  to  generate  the  snapshot  ensemble. 

As  before,  we  now  project  the  full  order  solution  at  each  time  step  onto  the  lumped  and  split  POD  bases. 
The  resulting  temporal  coefficients  are  compared  to  those  predicted  by  the  reduced  order  models.  The 
results  for  the  lumped  method  are  shown  in  Figure  8.  Split  method  results  are  shown  in  Figure  9.  As  seen  in 
those  figures,  the  accuracy  of  the  split  method  is  better.  For  temporal  coefficients  with  significant  frequency 
content,  values  predicted  by  the  reduced  model  are  virtually  identical  to  those  obtained  by  projecting  the  full 
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Lumped 


Split 


Figure  7:  Boundary  Condition  Accuracy. 


solution  onto  the  split  basis,  even  though  the  boundary  inputs  specified  are  different  than  those  used  during 
ensemble  creation.  The  split  method  is  better  suited  with  regard  to  feedback  control  law  design  as  it  is  more 
capable  of  accurately  representing  dynamics  that  are  not  explicitly  included  in  the  snapshot  ensemble. 


Lumped 


Figure  8:  Temporal  Coefficient  Accuracy  for  the  Lumped  Method. 


Closed-Loop  Results 

We  now  compare  the  effectiveness  of  the  lumped  and  split-POD  basis  methods  in  a  feedback  control 
setting.  The  reduced  system  given  by  (25),  (26)  is  linearized,  yielding  a  state-space  equation  of  the  form 

a(t)  =  Aa  +  Bu ,  (35) 

a(0)  =  «o-  (36) 

We  consider  the  tracking  problem  for  (35)-(36).  A  fixed  reference  signal  wref(x )  is  specified  for  the 
full  order  system.  Projecting  uye/(x)  onto  the  POD  basis  yields  tracking  coefficients  for  the  reduced  order 
model,  denoted  by  aref. 

To  formulate  the  tracking  control  problem,  we  consider  the  7-shifted  linear  quadratic  regulator  (LQR) 
cost  function 

/»oo 

J(a0,u)  =  /  {(a  —  aref)TQ(a  —  aref)  +  uT Ru]  e2ltdt.  (37) 

J  0 
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Split 


Figure  9:  Temporal  Coefficient  Accuracy  for  the  Split  Method. 


In  (37),  Q  is  a  diagonal,  symmetric,  positive  semi-definite  matrix  of  state  weights.  R  is  a  diagonal,  symmetric, 
positive  definite  matrix  of  control  weights.  The  quantity  7,  a  nonnegative  constant,  is  an  additional  parameter 
that  provides  added  robustness  in  the  control.27,28 

We  use  the  LQR  formulation  in  (37)  to  compare  closed-loop  results  obtained  from  the  lumped  and  split- 
POD  basis  methods.  A  snapshot  ensemble  is  constructed  containing  baseline  solution  data  as  well  as  data 
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Figure  10:  Split  Basis  Modes. 

resulting  from  nonzero  boundary  actuation.  As  it  is  desired  that  the  POD  basis  spans  unknown  dynamics 
introduced  by  the  LQR  feedback  control,  boundary  inputs  specified  during  ensemble  creation  are  of  the  form 

us{t)  =  Csin  ^7 r  (e4)0'3^  ur(t)  =  0,  (38) 

unit)  =  0  uT(f)  =  Csin  ^7r  (e4)°'3^  (39) 


for  C=  1,2,3. 

For  each  control  input  listed  in  (38)-(39),  snapshots  are  taken  in  increments  of  At  =  0.1  starting  from 
t  =  0  and  ending  at  t  =  15.  The  steady  baseline  solution  is  specified  for  the  initial  condition.  The  resulting 
snapshot  ensemble  consists  of  roughly  900  snapshots.  Requiring  that  99%  of  the  ensemble  energy  be  contained 
in  the  POD  basis  results  in  a  lumped  basis  consisting  of  5  modes.  A  split  basis  comprised  of  1  baseline  mode 
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and  20  actuator  modes  is  needed  when  baseline  and  actuated  energy  are  considered  separately.  The  first  9 
modes  of  the  split  basis  are  shown  in  Figure  10.  In  that  figure,  mode  1  is  the  baseline  mode.  Modes  2-9  are 
actuator  modes. 

The  tracking  LQR  problem  requires  the  specification  of  the  reference  signal  aref.  In  the  results  that 
follow,  aref  is  obtained  from  the  unactuated  steady  solution  for  the  case  Re  =  50.  This  solution  is  projected 


Figure  11:  Tracking  Reference  Function. 

onto  the  lumped  and  split  bases.  The  temporal  values  obtained  are  used  as  tracking  coefficients  in  the  reduced 
order  control  problem.  The  reference  signal  obtained  by  projecting  the  steady  solution  at  Re  =  50  onto  the 
split  basis  is  shown  in  Figure  11.  The  reference  function  obtained  by  projecting  onto  the  lumped  basis  is 
similar.  To  complete  the  control  formulation,  each  state  in  the  reduced  order  model  is  prescribed  a  weight 
of  2,500.  The  two  boundary  controls  are  each  given  unit  weight.  The  value  specified  for  7  in  (37)  is  0.25. 
The  closed-loop  solution  of  the  reduced  order  model  constructed  with  the  split  POD  basis  is  shown  in  Figure 
12.  By  comparing  the  controlled  solution  of  Figure  12  to  the  reference  function  in  Figure  11,  it  is  apparent 


t  =  0.05 


t  =  2.55 


Figure  12:  Closed-Loop  Split  Model. 

that  the  closed-loop  reduced  order  model  satisfies  the  control  objective  quite  well.  Separate  consideration  of 
actuated  energy  in  the  split  POD  basis  method  results  in  satisfactory  tracking  of  the  reference  signal. 

When  the  energy  ratio  in  (3)  is  applied  to  baseline  and  actuated  data  lumped  together  into  an  overall 
snapshot  ensemble,  the  results  are  much  less  favorable.  Closed-loop  solutions  of  the  reduced  order  model 
constructed  with  the  lumped  POD  basis  are  shown  in  Figure  13.  As  seen  in  that  figure,  virtually  no  tracking 
is  achieved  by  the  reduced  order  control.  Adjusting  parameters  in  the  control  formulation  has  little  effect 
on  this  result.  Increasing  the  state-weights  and  the  parameter  7  to  10,000  and  0.75,  respectively,  does  not 
significantly  improve  the  performance  of  the  control.  System  information  relevant  from  a  control  standpoint 


12 


t  =  0.05 


t  =  2.55 


Figure  13:  Closed-Loop  Lumped  Model. 


is  discarded  when  an  energy  argument  is  applied  during  order  reduction  to  the  lumped  snapshot  ensemble 
containing  baseline  and  actuated  data.  The  resulting  reduced  order  model  does  not  adequately  describe 
dynamics  induced  by  the  control.  Consequently  designing  a  feedback  control  from  such  a  model  results  in 
very  ineffective  response  when  the  control  is  applied  to  the  system. 

Full  Order  Validation 


To  validate  the  effectiveness  of  the  reduced  order  control  obtained  via  the  split  basis  method,  we  utilize  a 
fixed-point  projection  algorithm16  to  incorporate  the  reduced  order  boundary  control  in  the  full  order  model. 
The  closed-loop  solution  of  the  full  order  system  is  shown  in  Figure  14.  As  seen  in  that  figure,  the  reduced 
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Figure  14:  Closed-Loop  Response  with  Split  Method  Feedback  Control. 

order  control  effectively  drives  the  full  order  plant  to  the  target  profile.  The  full  order  discretized  model 
is  comprised  of  roughly  2,000  states.  The  reduced  model  obtained  via  the  split  basis  method  describes  21 
states.  As  a  result,  system  dimension  is  reduced  by  roughly  two  orders  of  magnitude  with  the  resulting 
reduced  order  control  being  quite  effective. 

Conclusions 


In  this  paper,  a  reduced  basis  construction  method  was  developed  allowing  for  separate  consideration 
of  baseline  and  actuated  dynamics  in  the  reduced  modeling  process.  A  prototype  initial  boundary  value 
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problem,  governed  by  the  two-dimensional  Burgers  equation,  was  formulated  to  demonstrate  the  utility  of 
the  method.  When  actuated  energy  was  considered  separately,  much  better  agreement  was  seen  between 
open-loop  solutions  of  the  reduced  and  full  order  systems.  Separate  consideration  of  energy  induced  by 
the  boundary  control  resulted  in  effective  feedback  control  for  the  reduced  and  full  order  systems.  When 
actuated  energy  was  not  explicitly  accounted  for  in  the  reduced  modeling  process,  the  resulting  feedback 
control  was  completely  ineffective  when  applied  to  the  system. 

These  results  demonstrate  the  need  for  separate  consideration  of  baseline  and  actuated  energy  in  the 
reduced  modeling  process  when  the  resulting  model  is  to  be  used  for  feedback  control  law  design.  Basis 
construction  relying  on  an  energy  argument  applied  to  a  lumped  snapshot  ensemble  containing  baseline  and 
actuated  data  can  result  in  important  control  information  being  discarded  during  order  reduction.  This  is 
particularly  the  case  in  a  boundary  control  setting  where  ensemble  energy  is  typically  dominated  by  that  in 
the  baseline  data.  Feedback  controls  developed  from  the  resulting  reduced  model  are  likely  ineffective  when 
applied  to  the  system.  Separate  consideration  of  dynamics  induced  by  boundary  control  input  results  in 
reduced  order  controllers  that  are  much  more  effective  when  applied  to  the  reduced  and  full  order  systems. 

In  a  companion  paper,20  we  apply  these  methods  in  an  experimental  laser  turret  application  with  the 
control  objective  of  regularizing  the  unsteady  flow  over  the  turret. 
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